HST Time-Series Photometry of the Transiting 
Planet of HD 209458 1 



Timothy M. Brown 2 , David Charbonneau 3,2 , Ronald L. Gilliland 4 , Robert W. Noyes 3 , 

and 

Adam Burrows 5 
timbrown@hao . ucar . edu 
dcharbonneauOcf a. harvard. edu 
rnoyes@cf a. harvard. edu 

gillilOstsci . edu 
aburrows@as . arizona . edu 

ABSTRACT 

We have observed 4 transits of the planet of HD 209458 using the STIS spectrograph 
on HST. Summing the recorded counts over wavelength between 582 nm and 638 nm 
yields a photometric time series with 80 s time sampling and relative precision of about 
1.1 x 10~ 4 per sample. The folded light curve can be fit within observational errors 
using a model consisting of an opaque circular planet transiting a limb-darkened stellar 
disk. In this way we estimate the planetary radius R p = 1.347 ± 0.060 Rj U p, the 
orbital inclination i = 86.68° ± 0.14°, the stellar radius R* = 1.146 ± 0.050 R e , and 
one parameter describing the stellar limb darkening. Our estimated radius is smaller 
than those from earlier studies, but is consistent within measurement errors, and is 
also consistent with theoretical estimates of the radii of irradiated Jupiter-like planets. 
Satellites or rings orbiting the planet would, if large enough, be apparent from distor- 
tions of the light curve or from irregularities in the transit timings. We find no evidence 
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for either satellites or rings, with upper limits on satellite radius and mass of 1.2 R® and 
3 Mq), respectively. Opaque rings, if present, must be smaller than 1.8 planetary radii 
in radial extent. The high level of photometric precision attained in this experiment 
confirms the feasibility of photometric detection of Earth-sized planets circling Sun-like 
stars. 

Subject headings: binaries: eclipsing - planetary systems - stars: individual (HD 209458) 
- techniques: photometric 

1. Introduction 

The low-mass companion to HD 209458 is the first extrasolar planet found to transit the disk 
of its parent star (Charbonneau, Brown, Latham, & Mayor 2000; Henry, Marcy, Butler, & Vogt 
2000). The primary star (GO V, V = 7.64, B-V = 0.58; H0g et al. (2000)) lies at at distance of 47 
pc as determined by Hipparcos (Perryman et al. 1997). An analysis of radial velocity measurements 
by Mazeh et al. (2000) gave an orbital period of 3.524 days, with M p sm.i = 0.69 ± 0.05 Mj up and 
a = 0.0468 AU, using the derived value of l.lztO.l M & for the stellar mass. When combined with the 
early photometric light curve data, the same analysis yielded an orbital inclination i = 86.1° ±1.6°, 
and a planetary radius R p = 1.40 ±0.17 -Rj U p- The planetary radius is at once the most interesting 
and the most uncertain of these parameters, largely because of uncertainty in the value of the 
stellar radius i?*. Knowledge of R p is important because it allows inferences about the planet's 
composition and evolutionary history (Guillot et al. 1996; Guillot 1999; Burrows et al. 2000). 
Unfortunately, the measured quantity that emerges most easily from the photometric transit data 
is the ratio R p /R*, and residual errors in the astrometry and effective stellar temperature suffice to 
make the estimate of R*, hence R p , uncertain by about 10%. Additional small errors in R p result 
from uncertainties about the stellar limb darkening. Jha et al. (2000) used multicolor photometric 
data to reduce these uncertainties, obtaining R p = 1.55±0.10 -Rj up - Finally, analyses of Hipparcos 
photometric data by Castellano et al. (2000) and Robichon & Arenou (2000) have given refined 
estimates of the orbital period. 

Here we report the results of very precise photometric measurements of transits of HD 209458 b, 
obtained using the STIS spectrograph on HST. The motivations for the study were 

1. to obtain sufficiently accurate photometry to reduce the ambiguity between estimates of the 
stellar radius, planetary radius, orbital inclination, and stellar limb darkening, 

2. to search for evidence of planetary satellites or circumplanetary rings, and 

3. to search the stellar spectrum observed in and out of transit for features imposed by trans- 
mission of starlight through the outer parts of the planet's atmosphere (Seager k, Sasselov 
2000; Brown 2001). 
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We discuss the photometric results (1 k, 2) in the present work; the spectroscopic investigation (3) 
will be the subject of a later paper. 



During a transit of HD 209458 b, the apparent brightness of the star is reduced by a little 
less than 2% as a result of the light blocked by the gas-giant planet. To have a useful sensitivity 
to smaller objects such as circumplanetary rings or Earth-sized satellites, we required photometric 
accuracies at least 2 orders of magnitude better, with time resolution of a few minutes or less. Most 
instruments on HST cannot meet these requirements, because they are not designed to accept the 
requisite large photon fluxes. The Space Telescope Imaging Spectrograph (STIS), used in spectro- 
graphic mode with the CCD detector, is the exception. By dispersing the light in a large bandwidth 
over many pixels, using the 4 e~ per data number gain setting, STIS can collect 2.5 x 10 8 photons 
per detector readout without saturating any detector pixels. A summation over all sampled wave- 
lengths then provides a photometric signal whose shot-noise-limited signal-to-noise ratio (SNR) is 
above 10 4 . The readout time for a detector subarray of the necessary size (1024 x 64 pixels) is 
about 20 s; it is thus possible to maintain both an 80 s sampling cadence and a respectable duty 
cycle of 75%. To attain near shot-noise limited operation, one must control, or verify the smallness 
of, numerous sources of instrumental noise, such as variations in the shutter opening and closing 
times, and CCD gain variations combined with spectrum motion on the detector. One of us had 
already investigated the use of the STIS as a precise photometer, and had demonstrated that, with 
care, these systematic errors can be suppressed (Gilliland 1999; Gilliland, Goudfrooij, & Kimble 
1999). We note that these studies further demonstrated that time series with SNR = 10 4 could be 
maintained even well past saturation of individual pixels, although in the present observations we 
did not saturate the detector. 

We observed HD 209458 during each of 4 planetary transits, on UT dates 2000 April 25, 
April 28/29, May 5/6, and May 12/13. The full duration of a transit is 184.25 minutes, or slightly 
less than 2 HST orbits, which are 96.5 minutes each. Objects orbiting the planet may, however, 
precede or follow it by as much as Rh, the radius of the planet's Hill sphere (i.e., the radius at 
which the star's tidal forces would dislodge a satellite): 



where a denotes the planet's orbital radius and M p and M* the masses of the planet and star. Since 
the planet's orbital speed is estimated to be 143 km s _1 , one must search for gravitationally-bound 
objects as much as 49 minutes, or about one half of an HST orbit, before and after the transit central 
time. This requires 4 full HST orbits to assure adequate coverage before and after the planetary 
transit. To allow the telescope pointing and thermal environment to stabilize before beginning 
critical observations, we added one further orbit at the beginning of each transit sequence, for a 



2. Observations and Data Analysis 




(1) 
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total of 5 HST orbits for each transit. The orbital phasing was such that each observed transit 
contained 2 initial orbits that were completely off-transit, 2 orbits that fell almost entirely during 
the transit, and one following orbit that was off-transit. 

We took all observations using the G750M grating, covering the wavelength range 581.3 < 
A < 638.2 nm, with a resolution of R = A/AA = 5540, corresponding to a resolution element of 
2.0 pixels or 0.11 nm. This wavelength range was chosen to cover the region of the Na D lines, where 
a large signature from transmission through the planetary atmosphere is thought to be likely. We 
used a 52 x 2 arcsec slit, much larger than the (typically 0.07 arcsec full width at half maximum) 
stellar image size, to minimize variations in the fraction of stellar light lost at the slit edges due to 
minor guiding and focus changes. (We did not use the "clear" slit in order to avoid excessive sky 
background.) The CCD area read out consisted of a 1024 x 64 subarray, covering the entire possible 
range in the dispersion direction and about 3.2 arcsec across the dispersion. Exposure times were 
60 s, which with a 20 s readout time gave an observing cadence of 1 sample per 80 s. This exposure 
time gave about 1.55 x 10 8 detected photons per spectrum, which would correspond to an optimal 
photon-noise-limited precision of 8.0 x 10~ 5 , or 87 /umag. Figure 1 shows a typical extracted stellar 
spectrum (see the discussion below for details of the extraction process). The first orbit of each 
group of 5 has a reduced number of spectra due to time spent on target acquisition. During each 
visit of five orbits, we obtained 28 spectra during the first visit, then 36 spectra during each of 
the subsequent three orbits, then 35 spectra during the final orbit, for a total of 171 spectra. The 
final data set thus contained 4 x 171 = 684 individual spectra of HD 209458. We took wavelength 
calibration spectra just before each of the 5 orbits, and, as an associated calibration program, we 
obtained 49 flat-field exposures during occultation time within the orbits surrounding the transit 
observations. All data for both the science program (8789) and calibration program (8797) are 
publically available via the HST archive 6 ; there is no proprietary period for these data since they 
were obtained through Director's Discretionary time. 

Observations of the first transit (UT 2000 April 25) were partly compromised by a database 
error in the location of the subarray for this rarely-used secondary central wavelength setting of 
G750M. The result was that the spectrum was not entirely contained within the CCD subarray 
for the red half of the spectral range. This error led to reduced counts and increased sensitivity 
to cross-dispersion positional drifts over this range. In the subsequent data analysis, we ignored 
this part of the wavelength range, but only for the first transit. Rapid work by the staff at STScI 
identified the source of this error and provided an effective correction strategy (a real-time slew 
executed after target acquisition for the visit on April 28/28, and a correction to the database 
thereafter), so that the spectra for subsequent transits were correctly positioned on the subarray. 

Data reduction consisted of 
1. recalibrating the 2-dimensional CCD images, 



See http://archive.stsci.edu. 
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2. removing cosmic ray events, 

3. extracting 1-dimensional spectra, 

4. summing the detected counts over wavelength to yield a photometric index 

5. correcting the resulting photometric time series for variations that depend on the phase of 
the HST orbit, and for variations between visits 

Recalibration of the CCD images started from the bias- and dark-subtracted and flat-fielded 
images produced by the standard STIS data reduction pipeline. The STIS pipeline builds a weekly 
bias image from exposures taken for this purpose. We examined the biases that were subtracted 
from each image and found them to be sufficiently accurate. In contrast to this, the pixel-to-pixel 
flat field images are often produced from data that are several months old, and we found that 
the relative sensitivity of numerous pixels had changed since the most recent pipeline flats were 
produced. Furthermore, the pipeline flats did not have a sufficiently high SNR for our purposes. 
We requested that high SNR flats be taken during times of Earth-occultation contemporaneous 
with our science exposures. We shot two sets of 7 such exposures for each of the visits on UT 2000 
April 24, 27, & May 12, and one set of 7 such exposures for the visit on UT 2000 May 5. For 
each set of 7 exposures, the intensity in each exposure increased with time and in a reproducible 
manner from one set to the next. The origin of this increase is unknown, although a possibility 
is that the calibration lamp warmed slowly as the series of exposures progressed. This effect was 
removed by renormalizing by the total number of counts. The relative illumination among pixels 
by the continuum source was remarkably constant: No temporal drift was observed (aside from 
the normalization variation described above), thus we were able to quantify the relative sensitivity 
of individual pixels (which typically varied with an amplitude of 1%) to a precision of better than 
than 0.1%. Between each visit, a small fraction of the pixels changed their relative sensitivity 
significantly, most likely due to radiation events between visits. Furthermore, the STIS CCD was 
annealed on UT 2000 May 6, between the third and fourth visits, which changed the behavior of 
a number of pixels. To handle sensitivity variations between any two visits, we tagged all pixels 
whose sensitivity we observed (from the contemporaneous flats) to have changed, and produced 
new flat fields for each visit where these pixels were replaced by their value derived from the data 
obtained during only the one visit. We then multiplied all the science exposures by the pipeline 
flat field, to undo the calibration from the pipeline, and then divided all the exposures by our 
newly derived flat fields. The calibration changes we imposed were typically small, typically a 
fraction of a percent in the gain of individual pixels. Since the photometric signal (described 
below) is produced by summing over a large number of pixels, these changes made only a small, 
but noticeable, improvement in the final results. In contrast to this, the improved precision from 
this recalibration provided significant gains in the analysis of the spectra to search for transmission 
features (to be presented in a separate paper). 

Cosmic ray hits on the CCD were a significant source of noise in the recalibrated data. We 
corrected for these by forming a time series of the intensity at each pixel and fitting this series as 
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a fourth-order polynomial in time. Points showing statistically significant positive differences from 
the fit were deemed to be corrupted. Such points were replaced by the value of the fitted function 
at the corresponding time. Of the 65536 pixels in each image, typically 36 (this is the median 
number) had to be corrected in this way, and of these pixels about 25% contributed to the final 
photometric index for the image. 

We extracted 1-dimensional spectra simply by summing along CCD columns, taking at each 
column a band of fixed width centered on the measured cross-dispersion spectrum position. Since 
motion of the spectrum in the cross-dispersion direction was measured to be much less than one 
pixel (typically motions were roughly 0.05 pixels within visits) for the entire data set, we used the 
same position for the integration band to produce all extracted spectra. We used a total cross- 
dispersion band 17 pixels in width, which we found produced the minimum rms variation in the 
photometric time-series for visits that occurred out of transit. The result of this operation was a 
1-dimensional spectrum sampled at 1024 points. 

For the most part we describe here a single photometric index, namely a sum S(t) of the 1- 
dimensional spectrum over almost the entire available spectral range (581.9 to 614.6 nm for the first 
transit, and 581.9 to 637.6 nm for the others). The only refinement in this process was to position 
the endpoints of the summation so as to avoid obvious absorption lines, so that the sum would 
be insensitive to displacements along the dispersion. In addition we performed an experiment in 
which we summed separately over the red and blue halves of the wavelength range; the difference 
between these two light curves results from the color dependence of the stellar limb darkening. The 
details of this analysis are described in section 3.3. 

The time series S(t) shows small but repeatable variations in phase with the HST orbital 
period. We do not know the origin of these variations, though we suspect that they are connected 
with the telescope's orbital thermal cycle. The variations are fairly well approximated by a linear 
decrease with an amplitude of 0.1 % over 48 minutes, with some curvature at the beginning and 
end of each orbit. We corrected the time series on a transit-by-transit basis. For each transit, 
we phased the data from the out-of-transit orbits (orbits 2 and 5) to the HST orbital period. We 
then fit a fourth order polynomial to these data, and divided all five orbits of the transit by this 
function. The first sample in every orbit is always smaller than the average by 0.25%; we rejected 
these values from the time series. 

The corrected time series have transit-to-transit differences of scale that are large for the first 
transit (because we summed over a smaller range of wavelength) and smaller (typically 0.1%) for 
the remaining 3 transits. Lacking a comparison star or other external calibration source, we cannot 
say how much of this long-term variation arises from the star and how much from instrumental drift. 
For our purposes this is not important, however, since we are primarily concerned with variations 
on time scales of minutes or hours, and not weeks. These scale variations were also removed by the 
procedure described in the previous paragraph. 

The normalized time series have a value of unity when averaged over the out-of-transit orbits, 
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and minimum values near the transit centers of 0.9836. Figure 2 shows the time series for each of 
the 4 transits phased with respect to the planetary orbit, and Figure 3 shows the combined data 
for all 4 phased transits, excluding the first orbit of each visit. One may estimate the noise in 
Figures 2 and 3 from the scatter among the out-of-transit observations. For the last 3 transits, 
these are typically 1.1 x 10~ 4 , or 120 /xmag for each 60 s integration; the noise for the first transit 
is worse by a factor of about 2.3, mostly because of increased sensitivity to transverse motions of 
the poorly-centered spectra. This precision suffices to show the 1.64% transit dip with a SNR of 
about 150, at high time resolution. 

3. Interpretation 

3.1. Planetary Orbital Parameters 

Ideally, we would like to measure T c , the time at the center of the transit, for each of the 
four observed transit events. However, we have gaps in our observations of any one transit (see 
Figure 2) due to Earth occultation that prevent us from doing so with sufficient accuracy. We 
can, however, estimate the period accurately (and independently from the absolute times of the 
transits) by seeing which assumed value of the period produces the minimum scatter in the phased 
transit curve. In this manner we derive an orbital period of P = 3.52474 ± 0.00007 d. Having 
evaluated the period, we can then subtract the appropriate multiple of it from each of visits 2, 3, 
& 4 (1, 3, & 5 planetary orbital periods later, respectively), to phase the data to the time of the 
first transit. We then fold the transit curve about an assumed time T c for the midpoint, and derive 
T c from minimizing the scatter between the observation before and after the center of transit. We 
derive T c = 2451659.93675 ± 0.0001 HJD. We can then compare this value of T c with the three 
other accurately measured transit times, two from Charbonneau, Brown, Latham, & Mayor (2000) 
and one from Jha et al. (2000). Doing so, we derive a more accurate value of the period, due 
to the large number of transits that have now elapsed from these observation from fall 1999. We 
find P = 3.52480 ± 0.00004 d. The uncertainty in the period is 4.3 seconds. Our derived value 
is in excellent agreement with the values derived with similar precision from Hipparcos archive 
photometry by Castellano et al. (2000), who found P = 3.524736 ± 0.000045 d, and by Robichon 
& Arenou (2000), who found P = 3.524739 ± 0.000014 d. 

3.2. Stellar and Planetary Parameters 

To estimate the planet's radius, we sought to represent the light curve as the result of an 
opaque sphere of radius R p in an inclined circular orbit about a limb-darkened star of radius R*. 
We took the stellar mass to be 1.1 ± 0.1 M Q (Mazeh et al. 2000), with the orbital radius and 
velocity determined from the period and Kepler's laws. Free parameters in the fit were R p , i?*, 
the orbital inclination i, and the parameters u\ and U2 describing the stellar limb darkening in the 
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form 

W =1 _ tii(1 _ Ai) _ U2(1 _ Ai) 2 j (2) 

where n is the cosine of the angle between the line of sight and the normal to the local stellar surface 
(Claret h Gimenez 1990). A simpler and often-used parameterization replaces the right hand side 
of Eq. (2) with 1 — u(l — /x). We find that this parameterization worsens the x 2 statistic of the fit 
slightly, but significantly. We note, however, that fitting the expression in Eq. (2) gives strongly 
anticorrelated formal errors for u\ and U2- Thus, a more natural set of limb darkening parameters is 
(ui + U2) and (tii — U2)- The former describes the total magnitude of the limb darkening variation, 
and is well constrained by the observations, while the latter describes the amount of curvature in 
the limb darkening function, and is relatively poorly constrained. For explanatory purposes below, 
we shall refer to these combinations rather than to u\ and 112 individually, and moreover we shall 
adopt the approximation that u ~ u\ + 112 ■ 

The conceptual basis for this fitting process is illustrated in Figure 4. At the lowest level of 
approximation, the light curve is described by just 2 parameters: its depth d depends mostly upon 
Rp/R*, while its duration I depends mostly upon the transit's chord length, and therefore upon i?* 
and i. At this treatment level there are fewer observables {d, 1} than unknowns {Rp, R*,i}, and 
so one must estimate R* from other evidence in order to obtain values for R p and i. This fitting 
degeneracy can be removed by taking account of more subtle effects: The duration w of the planet's 
ingress and egress depends upon R p but is also proportional to secV>, where ip is the angle between 
the planet's line of motion and the local normal to the stellar limb. Thus, w depends upon Rp, i?*, 
and i. Finally, the curvature C of the light curve between second and third contacts depends upon 
the stellar limb darkening parameter u and upon i and R*. Thus, if d, I, w, and C can be measured 
with adequate precision, one may estimate each of the 4 independent system parameters R*, Rp, 
i, and u. Fitting for both u\ and U2 requires, in addition to the quantities already mentioned, a 
measurement of the detailed shape of the light curve between second and third contacts. All of the 
foregoing assumes that the star's mass M#, and hence the planet's orbital velocity and semi-major 
axis, are known. The derived value of R p is only weakly dependent upon the assumed value of M*, 
scaling as R p cx M* 1 / 3 . We assume M* = 1.1 ± 0.1 M (Mazeh et al. 2000), and our derived 
errors include this uncertainty in the stellar mass. 

We derive best-fit values for {R p , R*, i, ui, U2} by minimizing the x 2 of the fit. The reduced 
X 2 for the best fit values was 1.07, indicating that the model is a good fit to these data. To derive 
la errors for each parameter, we change the value of that parameter and fix it at a new value, and 
then allow all other parameters to float, as well as allow for a stellar mass between 1.0 and 1.2 Mq. 
We repeat this procedure until the best fit solution produces an increase in the x 2 corresponding 
to a la change. The best-fit parameters and their errors are given in Table 1. Both the stellar 
and planetary radii are found to be in agreement with those derived by Mazeh et al. (2000), 
and the precision has been increased greatly. Most significantly, the error in the planetary radius 
has been reduced from 0.17 i?j up to 0.06 i?j up . We note also that our derived value for R* is in 
agreement with that derived from the Hipparcos distance, i?* = 1.18 ± 0.09i? Q (see Jha et al. 2000 
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for details). The values for R p and R* derived here are somewhat smaller than those derived by 
analysis of multicolor data by Jha et al. (2000), although the differences are consistent with the 
errors, differing by less than 1.4a. There are some correlations among the formal errors for the 
various derived parameters, and these are illustrated in Figure 5. We reiterate that our stated la 
errors include, and are dominated by, these correlations. In particular, the uncertainties in the 
radii of the planet and star are dominated by the assumed uncertainty of ±0.1 M Q in the stellar 
mass. 



3.3. Color Dependence of the Light Curve 

To check for color dependence of the transit curve, we divided the spectral region in two ranges, 
581.9-598.3 nm ("blue"), and 598.3-637.6 nm ("red"). We then generated photometric timeseries 
for each of these following the procedure described in Section 2. We used data from only visits 2, 
3, and 4, since the data from visit 1 lacked the red half. For each color, we phased the data to 
the period of 3.52474 d, and folded it about the midpoint of the transit, T c . Since we do not have 
useful data from the first visit, gaps exist in the time coverage. We then grouped the data into 
5-minute bins, and for the data in each bin, we fit brightness as a linear function of time within the 
bin. We took the brightness for the bin to be this linear function, evaluated at the centeral time 
for the bin. We then differenced the red and blue transit curves generated in this manner. These 
data are shown in Figure 6. From these data, it is clear that the transit is deeper in the blue at 
times near the center of transit, and deeper in the red at times when the planet is near the limb, 
as would be expected from the greater limb-darkening in the blue (see, e.g., Rosenblatt (1971); 
Sackett (1999) for sample differenced-color light curves). 

To generate the best-fit light curve for the difference of these color data, we fixed Rp, R*, and 
i at the values derived in Section 3.1. We then allowed small changes to the parameters describing 
the limb darkening, such that for the red half (denoted by R), 

(ui + u 2 ) R = (ui + u 2 ) + a 

(m - u 2 ) R = (u! - u 2 ) + (3, (3) 
and for the blue half (denoted by B), 

(ui + u 2 ) B = (ui +u 2 ) - a 

(m -u 2 ) B = (ui - u 2 ) - (3. (4) 

We then generated model light curves for both the red and blue data, differenced these, and 
evaluated the \ 2 of this fit to the data. The data are best fit by a = —0.021 and (3 = 0.063, thus 
the derived values for the limb-darkening are 

{u\+u 2 ) R = 0.619 ±0.03 
(«i + u 2 ) B = 0.661 ±0.03 
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(tii -u 2 ) R = 0.003 ±0.1 
(«i -u 2 ) B = -0.123 ±0.1, (5) 

consistent with the measured Solar limb-darkening (Cox 1999). We show this best-fit curve in 
Figure 6. Assuming the best-fit values for i and R*, we can calculate the projected separation of 
the center of the planet and the center of the star at any time during the transit, and thus at any 
point in this color curve. This axis is shown at the top of Figure 6. From this one can determine the 
fractional stellar radius at which the limb-darkening curves (weighted over the stellar disk) cross; 
this occurs near a radius of 0.84 R*. 

The color-dependent amplitude shown in Figure 6 is small primarily because we are contrasting 
the limb-darkening over two bands that are separated by only 28 nm. Observing the transit in very 
disparate bands across the visible and near-IR would show a much greater effect. If one were 
willing to assume a model for the limb-darkening of the star, one could use this effect to break 
the degeneracy between the parameters R p , i?*, and i, even if the data did not have the precision 
of that which we present here. This is precisely what was done with multicolor observations of 
HD 209458 in Jha et al. (2000). 



3.4. Search for Circumplanetary Rings 

If the planet of HD 209458 were circled by a ring system with significant opacity, the rings 
would cause distortions of the light curve relative to that of a spherical body (Schneider 1999). 
The cross-sectional area of the planet would appear larger because of the light obstructed by the 
rings, and (more usefully, for detection purposes) one would also see small dips in the light curve 
before first and after fourth contact (see Figure 7). The phased light curve is fit within observational 
errors by the simple planetary-transit model described above, so there is no evidence for rings in the 
current data. We can, however, set an upper limit on the size of a ring system consistent with the 
observations. For this purpose we assume that such a ring system lies in the planet's orbital plane, 
that it extends continuously from the planet's limb to a maximum radius R r (measured in units 
of Rp), and that it is entirely opaque to transmitted light. We cannot simply assume the best-fit 
values for {Rp, i?*, i, ui, u 2 } from section 3.1, because, if circumplanetary rings are present, some 
area is occulted by the rings; this would cause us to over-estimate R p . To simplify the investigation, 
we fix i?*, u\, and u 2 . Changing the value of -R* does not affect the results since it would result 
in a larger value of R p , and we state an upper limit for R r in units of R p . For each trial value of 
Rr, we allow R p and i to float, and derive an upper limit for R r by finding the value above which 
the x 2 increases by an amount corresponding to a 3a change. The maximum ring radius consistent 
with the observations at this confidence level is then 1.8 R p ; this is slightly smaller than the radius 
of Saturn's ring system, measured in units of Saturn's radius. This fairly low sensitivity to ring 
systems results mostly from the assumption that rings must lie in the planet's orbital plane, and 
hence that the ring plane must be nearly edge-on as seen from Earth. 
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3.5. Search for Planetary Satellites 

A satellite orbiting HD 209458 b might be detectable either from its photometric signature 
or from its influence on the orbital motion of the planet itself (Sartoretti & Schneider 1999). A 
satellite would block light in addition to that obstructed by the planet, unless it happened to be 
projected onto the planet's disk during the transit. This additional obstruction could occur either 
earlier or later than the main transit, depending upon the satellite's position in its orbit. Similarly, 
its duration could be slightly different from that of the planetary transit, because the satellite can 
move significantly in its orbit during the transit, and because the chord the satellite strikes across 
the star may be longer or shorter than that of the planet. 

As was the case for rings, there is no evidence for satellites in the photometric time series. 
Figure 8 shows the residuals about the best-fit light curve for each of the observed transits, along 
with transit curves that might be expected from a satellite with 1.5 times the radius of the Earth 
and an orbital period of 1.5 d. A signal of this size would be easily detected, if it were present. 
To set a better limit on the size of possible satellites, we searched the residual data for repeated 
transit-like events by applying matched filters that simulated the light curves from satellites with a 
range of orbital periods, phases, and semi-major axes. This process was analogous to that used to 
search for transiting planets in observations of the eclipsing binary CM Dra (Doyle et al. 2000), 
and in HST time-series photometry of the globular cluster 47 Tucanae (Gilliland et al. 2000; 
Brown et al. 2001). The process we actually used worked as follows: We first subtracted the 
best-fit model for the transit curve of the planet from the data, yielding the photometric residuals 
R(i), with i = [1, N]. We used the data only from the last three visits, as the errors due to 
the offset of the spectrum in the CCD subarray during the first transit are relatively very large. 
We constructed a dense grid of sample satellite orbital periods (p € [1 d, 3.5 d]) and satellite 
orbital phases {<p <G [0,1]). For each pair of {p,4>}, we evaluated, for each visit, the amount of 
time by which the center of the satellite transit would lead or trail the center of the planetary 
transit. We then subtracted these corrections from the times of the observations (so as to phase 
the data to any potential satellite transit), yielding the new times t P:< p(i). Neglecting small changes 
in the duration of the satellite transit due to satellite orbital motion during the transit times, we 
assumed the duration of the satellite transit, d sa t, to be the same as that of the planetary transit, 
184.25 minutes. We then computed the correlation of these data with a box-car function defined 
by: 

1 if |t p ^(t)| <d sat /2 1 
-1 if \t p j(i)\ > d sat /2 j ' 



B p, *(«) 

The correlation was given by 



c p,4> = ' W 

where the leading factor of 2 accounts for the fact that any such satellite transit will have a mean 
of zero, since we had already subtracted the best-fit planetary transit. In the case of noiseless data, 
the value of C p ^ at the correct {p, <f>} is the depth of the transit: For a satellite of radius R sa t 
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passing in front of HD 209458, we would have C p ^ ~ 6.4 x 10 b (R sat / R®) 2 ■ 

To derive detection thresholds, we replaced all the residuals with numbers drawn from a 
Gaussian distribution with a = 1.14 x 10~ 4 and evaluated C v ^. We repeated this many times, and 
evaluated the maximum value of C v ^ for each fake data set. These maximum values were Gaussian- 
distributed, with the upper 3-a point of the distribution lying at C = 5.2 x 10~ 5 . Thus, assuming 
the noise to be purely Gaussian, we exclude with 99.7% confidence the presence of satellites larger 
than the corresponding radius, namely R sa t = 0.9 JS®. 

The above conclusion is correct only if the errors are indeed Gaussian. However, we can see 
from the residuals that there may exist orbit-to-orbit drifts. To account for these as well, we 
performed a different test. For each orbit, we calculated the average of the residuals, and found the 
maximum average offset to be 6.9 x 10~ 5 . As before, we replaced all the residuals with numbers 
drawn from a Gaussian distribution with a = 1.14 x 10 -4 . We then added to each orbit of data 
an additional offset drawn from a uniform distribution between ±6.9 x 10~ 5 . We then evaluated 
C Pi( p, and repeated this procedure many times. Following the same procedure as before, we found 
a 3a limit of C = 9.0 x 10~ 5 , corresponding to R sat = 1.2 i?®. We believe that this noise model 
better describes our data than does pure Gaussian noise, so we prefer this estimate of the detection 
threshold. The smaller threshold remains of interest, however, as an indication of the size of 
detectable companions, should the remaining systematic errors of measurement be eliminated. 

For our actual residuals, the maximum of C p ^ is 4.88 x 10~ 5 . This value is significantly below 
the detection threshold described in the last paragraph, and consistent with the result of pure 
noise. We conclude that we have no evidence for a satellite orbiting the planet of HD 209458, and 
can exclude (at the 3a level) satellites larger in radius than 1.2 i?®. These results show that if we 
obtained the same precision on a star with i?* = 1.0 R Q (recall that for HD 209458, i?* = 1.15 Rq), 
we would be able to make a 3a detection of a 1.0 — -R®-sized planet in transit across its star. In 
the absence of orbit-to-orbit drifts, the limit would be 0.78 i?®. 

Satellites might also be detected by the periodic displacement of the planet in its orbit due 
to the gravitational attraction of the satellite. The magnitude 6x of this displacement could be as 
much as Sx = a sa tM sa t/M p , where a sa t is the satellite's orbital radius and M sa t is its mass. The 
visible effect would be to wobble the planet ahead of or behind its mean orbital phase, assuming the 
satellite orbit to be approximately coplanar with that of the planet. Transits would therefore occur 
early or late relative to the ephemeris, depending upon the phase of the satellite in its orbit. For 
an Earth- mass satellite orbiting HD 209458 b at a distance of one Hill sphere radius, the maximum 
temporal excursion is 13 s. This time is comparable to the formal 3a error on the estimate of the 
central time of a single transit. Relative to the time predicted from a best-fit estimate (based on 
these HST data only) of the period and initial epoch, the observed timing displacements (seconds) 
for the 4 transits are {-20.0 ±10.0, 0.3±4.4, -3.3±6.2, 10.0±4.6}. Figure 9 shows these observed 
displacements. Also shown is the % 2 statistic for each transit, as a function of displacement from the 
observed transit times. From these curves, it is evident that the differences between the observed 
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transit times and the predicted ones are at most about 2a. We doubt these discrepancies arise 
from satellites of the planet. Rather, we suspect that small systematic errors in the observations 
(for example, small linear trends in detector sensitivity that run from beginning to end of a 5-orbit 
transit observation) may be responsible. The residuals shown in Fig. 8 support this conclusion, 
with the residuals for the transit on 12 May being predominantly negative before the transit and 
positive after it; residuals due to a planetary satellite should affect the in-transit data (orbits 3 and 
4), but not the out-of-transit data (orbits 2 and 5). Taking the observed transit time variations 
to be an upper limit on the displacement caused by an unseen satellite, and assuming its orbital 
radius to be that of the Hill sphere, we can exclude with 3a confidence the presence of satellites of 
more than 3 Earth masses. 



4. Discussion 

The planetary radius 1.347 Rj up inferred from the HST light curve is consistent with a pre- 
viously reported value of 1.40 .Rj^ (Mazeh et al. 2000; Henry, Marcy, Butler, & Vogt 2000), 
and somewhat smaller than the value of 1.55 Rj up derived from multi-color data by Jha et al. 
(2000). The new estimate is more likely to be correct, because the radius of the parent star is 
determined as part of the fitting process, rather than being assumed. Indeed, the works just cited 
identified uncertainty in the stellar radius as an important contributor to the uncertainty of the 
final result. An assumption about the stellar mass is required in any case, but our errors include 
a 10% uncertainty in the stellar mass. Furthermore, the derived value for the planetary radius 
is only weakly dependent upon the value of the stellar mass. A radius of 1.35 i?j up is consistent 
with the irradiated model described by Burrows et al. (2000); it falls between the models with 
high (0.5) and low (0.0) Bond albedo, and the error of 0.05 Rj up is, in principle, small enough to 
distinguish between these two albedos. Implied values of the mean density p, surface gravity g, 
and escape velocity v e all increase relative to earlier estimates, but not by large amounts. We find 
p = 0.35 gem -3 , g = 943 cms -2 , and v e = 43 kms -1 . These changes are all such as to increase the 
estimated stability of the planet against disruption by tidal forces, thermal evaporation, or mass 
stripping by the stellar wind. 

Although it was possible a priori for the system to have observably large satellites or rings, 
the absence of these features is not surprising. Only large satellites (bigger than the Earth) could 
have been detected by these observations, and for a satellite to survive so close to the star, it 
would have to be made of refractory materials. The solar system contains no bodies that meet 
both of these requirements (Earth and Venus come the closest), so perhaps it is reasonable to 
guess that HD 209458 b likewise is not home to such an object. Similar comments apply to a ring 
system, which would have to be large and opaque in order be detectable in our observations. In 
the harsh radiation environment 0.05 AU from the central star, the processes that destroy rings of 
fine particles would be accelerated, and long lifetimes may not be expected. One must remember, 
however, that the present observations could not have detected the Galilean satellites of Jupiter, 
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and (because of their small optical depth) the rings of Saturn would be only marginally visible. 
Our observations constrain the presence of such objects only to the extent that large and obvious 
companions are excluded. 

In addition to their direct scientific interest, these observations provide the best example to 
date of the capabilities of extremely precise photometry from space. This suggest a number of 
potentially rewarding future observations: 

1. With the achieved precision, it may be possible to detect the reflected light from a close-in 
giant planet (Seager & Sasselov 1998; Sudarsky, Burrows & Pinto 2000; Marley et al. 1999). 
To date, such studies (Charbonneau et al. 1999; Collier Cameron et al. 1999) have produced 
only upper limits, and have required restrictive assumptions about the spectrum of the reflected 
light. The amplitude of this effect in the HD 209458 system would be e\ = p\(R p /a) 2 = p\ x 
2.0 x 10 -4 , where p\ is the wavelength-dependent geometric albedo. The transiting configuration 
of HD 209458 b would make it the ideal target, since at times of secondary eclipse (ie. just as the 
planet passes behind the star), the planet would pass within 26 minutes from being nearly fully 
illuminated to not visible. For such observations where there is no particular interest in moderate 
resolution spectroscopy (as was the case for this program), it would be desirable to switch to the 
low resolution gratings, and thus measure the albedo across the wavelength range where the central 
star outputs the majority of its energy. One could thus evaluate the net energy deposition into the 
planet, a key quantity in understanding its evolution (Burrows et al. 2000). 

2. The perturbation caused by additional planets in the HD 209458 system would change 
the observed times of transit. A 1 Mj up planet at 10 AU would cause the central star to move 
± 0.01 AU, and thus the transits would be observed as much as 5 s earlier or later. This effect 
could be used to infer the presence of such additional companions, although many years would be 
required to observe the effect, due to the long orbital periods at these large semi-major axes. 

3. Although the STIS instrument was not designed with high SNR photometry in mind, the 
achieved precision of the photometric time series confirms that it is feasible to detect the transits 
of Earth-sized planets across the disks of Sun-like stars. Moreover, reaching this precision does not 
depend upon the HST's large aperture; the limited bandwidth used in this experiment and (to a 
lesser degree) the transmission losses in the spectrograph combined to make the photometry far 
less efficient than it could be. For example, with a bandwidth spanning 400-1000 nm and a system 
efficiency of 50%, the photon count rates achieved here could be reached with a telescope of only 
25 cm aperture. Placed in a suitable orbit, so that full uninterrupted transits of HD 209458 b 
could be observed, such a telescope could detect a satellite of the mass or radius of Ganymede after 
observing about 100 transits, which is to say, within about 1 year. The first generation of orbiting 
telescopes designed for such purposes are now under development (MOST, Matthews et al. (2000); 
COROT, Michel et al. (2000); MONS, Kjeldsen, Bedding, & Christensen-Dalsgaard (2000)), and 
larger-scale projects are planned or proposed, such as The Kepler Mission, (Koch et al. 1998), and 
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Eddington 7 . 

With the new photometric and spectroscopic data sets that will soon become available, we can 
look forward to an exciting decade unraveling the structure and history of the close-in extrasolar 
planets. 

We are grateful to the HST STIS and operations teams, especially Helen Hart, Gerard Kriss, 
and Jeff Valenti for their prompt and insightful help in resolving the database error experienced dur- 
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a grant from the Space Telescope Science Institute, which is operated by the Association of Uni- 
versities for Research in Astronomy, Inc., under NASA contract NAS5-26555. Further support for 
this work was provided through NASA grants NAG5-7073 and NAG5-7499. D. C. is supported in 
part by a Newkirk Fellowship of the High Altitude Observatory. 
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Table 1. HD 209458 b Fit Parameters 



R* 1.146 ± 0.050 R { 

R p 1.347 ± 0.060 R 3 . 

i 86.68° ± 0.14° 

(m + u 2 ) 0.640 ± 0.030 

(u!-u 2 ) -0.055 ± 0.100 
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Fig. 1. — Typical spectrum of HD 209458 obtained with the STIS instrument, showing the number 
of detected photons in a single 60 s integration vs wavelength. We obtained 684 such spectra in all. 

Fig. 2. — Time series of the corrected intensity shown separately for each of the 4 transits observed 
by HST, with successive transits offset by -0.006 for clarity. Note that, because the transit duration 
is almost 2 HST orbits, complete temporal coverage was not obtained for any one transit. 

Fig. 3. — Phased light curve for all 4 transits, assuming a planetary orbital period of 3.52474 d. 
The time series for each transit has been scaled to have the same average intensity over the 2nd 
and 5th (out-of-transit) orbits. 

Fig. 4. — Schematic illustration of the light curve of a transiting planet. Measurable quantities are 
the duration of the transit I, the transit depth d, the ingress/egress duration w, and the central 
curvature of the light curve C. Given the orbital speed (which follows from the orbital period and 
the stellar mass), these quantities determine the radii of the star and of the planet, the orbital 
inclination, and the degree of limb darkening. 

Fig. 5. — The x 2 surfaces for the fit of R p , i?*, i, and (u\ + U2) to the observations. Each panel 
shows one 2-dimensional cut through the 5-dimensional x 2 surface, as indicated in the Figure. 
The displayed ranges in each plot are the la error intervals in each case. The plotted contours 
serve to indicate the mutual dependences of the selected parameters, but the contour levels do not 
correspond to specific significance levels. 

Fig. 6. — Color dependence of the transit light curve, shown as the difference (blue-red) between 
the normalized fluxes in the wavelength ranges [581.9 nm, 598.3 nm] (blue) and [598.3 nm, 637.6 
nm] (red), plotted against time from transit center. The solid line corresponds to the best-fit model 
color-dependent limb darkening law, as described by Eq. (5). 

Fig. 7. — Upper panel: the observed light curve folded around the central transit time (bars) binned 
to 60 s sampling. The bar height indicates the formal errors of the binned estimates. Regions 1 
and 2 are shown at larger scale in the lower panels. Lower panels: Observed intensities (black 
diamonds with la error bars) and curves for best-fit models assuming maximum ring radii of 1 R p 
(i.e., no rings) (solid curve), 1.8 R p (dotted curve) and 2.2 R p (dashed curve). The latter 2 curves 
are inconsistent with the data. Rings were assumed to be opaque, to lie in the planet's orbital 
plane, and to extend from 1 R p to the stated radii. 

Fig. 8. — Residuals around the fit of a 5-parameter transit model (see text) to the observations, 
shown with error bars. The solid lines show the light curve one would expect given the presence 
of a satellite of 1.5 Earth radii, with an orbital period of 1.5 d. Note the relatively large errors for 
the first transit, due to the offset of the spectrum in the CCD subarray. 

Fig. 9. — The \ 2 curves for fits to the central time of transit for each of the 4 transits (dotted 
= 25 April, dashed = 28/29 April, triple-dot dashed = 5/6 May, dot-dashed = 12/13 May). The 
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horizontal bar shows the range of timings resulting from a 3-Earth-mass satellite at a distance of 1 
Hill Sphere radius. 
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